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Abstract 

In this paper we present an algorithm for construction of minimal involutive 
polynomial bases which are Grobner bases of the special form. The most general 
involutive algorithms are based on the concept of involutive monomial division 
which leads to partition of variables into multiplicative and non- multiplicative. 
This partition gives thereby the self-consistent computational procedure for con- 
structing an involutive basis by performing non-multiplicative prolongations and 
multiplicative reductions. Every specific involutive division generates a particu- 
lar form of involutive computational procedure. In addition to three involutive 
divisions used by Thomas, Janet and Pommaret for analysis of partial differential 
equations we define two new ones. These two divisions, as well as Thomas divi- 
sion, do not depend on the order of variables. We prove noetherity, continuity 
and constructivity of the new divisions that provides correctness and termination 
of involutive algorithms for any finite set of input polynomials and any admissi- 
ble monomial ordering. We show that, given an admissible monomial ordering, a 
monic minimal involutive basis is uniquely defined and thereby can be considered 
as canonical much like the reduced Grobner basis. 



1 Introduction 

Computational aspects of constructing Grobner bases invented by Buchberger [[!]] are 
now under intensive investigation due to the great theoretical and practical importance 
of these bases in computational commutative algebra and algebraic geometry 0, |3|, [| . 
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Grobner bases are also becoming of greater importance in non-commutative U ^| 0] 
and differential algebra || [|. 

Since its invention about thirty years ago, feasibility of the Buchberger algorithm 
has been notably increased. First of all, it was resulted from discovering criteria 
for avoiding unnecessary reductions [|10| , |TT| , |T2"fl which allow a partial extension to 



non- commutative case 0. Next, the key role of the reduction and, especially, se- 
lection strategies was experimentally observed, and heuristically good strategies were 



found [13]. For construction of a lexicographical Grobner basis, which is the most use- 



ful for solving polynomial equations, an efficient computation scheme was developed 
in 



14] based on converting a basis from one ordering into another. 



On the other hand, Zharkov and Blinkov |15| were pioneered in revealing another 
computational scheme for Grobner bases construction in commutative algebra. They 
used the partition of variables into multiplicative and non-multiplicative invented in 
Pommaret (16| to bring partial differential equations into so-called involutive form [17] 
which has all the integrability conditions satisfied. Zharkov and Blinkov showed that 
sequential multiplication of the polynomials in the system by non-multiplicative vari- 
ables, and reduction of these prolonged polynomials modulo others, by means of their 
multiplicative power products only, ends up, under certain conditions, with a Grobner 
basis. Though the latter is generally not the reduced basis, it reveals some attractive 
features fl8| . 

Already first computer experiments carried out in [15j showed rather high efficiency 
of the new computational scheme. However, that algorithm terminates, generally, only 
for zero-dimensional ideals and for degree compatible term orderings JTP[ . The algebraic 
origin of such an algorithmic behavior was analyzed in 0] where it was also shown 
that Pommaret involutive bases are just Grobner ones of ideals in the commutative 
rings with respect to non-commutative gradings. Interconnection of Pommaret bases 
and Grobner bases was recently investigated also in [^] . 

In our previous paper |22] general algorithmic foundations of involutive approach to 
commutative algebra were considered, and a number of new concepts was introduced 
allowing one to study the involutive algorithmic procedure in its general form. The 
central concept of our analysis is involutive monomial division. Every specific invo- 
lutive division generates some particular computation procedure for constructing the 
corresponding involutive basis. Every involutive basis, if it is finite, was proved to be 
a Grobner basis, generally, redundant. We formulated the axiomatic properties of an 
involutive division which provide a proper partition of variables into multiplicative and 
non-multiplicative, and, hence, to construct different divisions. It was also proved that 
those partitions used by Janet [F7|, Thomas [23] and Pommaret |16j are generated by 
particular involutive divisions. 

Important properties of noetherity, continuity and construct ivity for an involutive 
division were also characterized. Noetherity provides for the existence of a finite invo- 
lutive monomial basis for any monomial ideal much like to the conventional monomial 
bases. Continuity assures involutivity of every locally involutive set. Constructivity is 
a strengthening of continuity. It allows one to compute an involutive monomial basis 
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from the initial one by means of its enlargement with single non-multiplicative pro- 
longations only, that is, to avoid enlargement with multiplicative prolongations. We 
showed that Janet and Thomas divisions are noetherian, continuous and constructive 
whereas Pommaret division, being continuous and constructive, is not noetherian. Just 
by this reason a positive-dimensional polynomial ideal, generally, does not have a finite 
Pommaret basis. We presented in |22] a general form of the involutive algorithm. Its 
correctness follows from continuity of a division while termination holds for any polyno- 
mial ideal and for any admissible monomial ordering only for noetherian divisions. The 
algorithm involves the Buchberger's chain criterion to avoid unnecessary reductions. 

In the present paper, in addition to Janet, Thomas and Pommaret divisions ana- 
lyzed in , we give examples of two more involutive divisions which are proved to be 
continuous, constructive and noetherian. We present also the special form of an invo- 
lutive algorithm which, given a constructive noetherian division, provides computation 
of a minimal involutive basis. We show that the monic form of the latter is uniquely 
defined for any fixed admissible monomial ordering. 

The rest of the paper is organized as follows. In Section 2 we give a brief review of 
involutive concepts and methods which are used in the following sections. In Section 3 
we consider some examples of involutive monomial divisions including those introduced 
by Thomas, Janet and Pommaret along with two new ones. In Section 4 we study 
the minimal involutive monomial bases. The algorithm for construction of minimal 
polynomial bases is described in Section 5, and some concluding remarks are given in 
Section 6. 



2 Background of Involutive Approach 

In this section we briefly describe the fundamentals of the general involutive approach 
proposed in [pi which are used in Sections 3-5. 

Let N be a set of non-negative integers, and M = {xf 1 ■ ■ ■ xf^ | dj G N} be a set of 
monomials in the polynomial ring R = K[x\, . . . , x n ] over zero characteristic field K. 

By deg(u) and degi(u) we denote the total degree of u G M and the degree of 
variable Xi in u, respectively. An admissible monomial ordering is denoted by and 
throughout this paper we shall assume that 

xi y x 2 y ■ ■ ■ y x n . (l) 

The leading monomial and the leading coefficient of polynomial / G R with respect 
to ordering -< are denoted by lm(f) and lc(f), respectively. If F C R is a polynomial 
set, then by lm(F) we denote the leading monomial set for F, and Id(F) will denote 
the ideal in R generated by F. For the least common multiple and for the greatest 
common divisor of two monomials u, v G M we shall use the conventional notations 
lcm(u,v) and gcd(u,v), respectively. 

If monomial u divides monomial v we shall write u\v. 



3 



Definition 2.1 An involutive division L on M is given, if for any finite monomial set 
U C M and for any u G U there is given a submonoid L(u,U) of M satisfying the 
conditions: 

(a) If w G L(u, U) and v\w, then t> G £7). 

(b) If u, v G £7 and uL(u, U) H uL(t>, C7) 7^ 0, then « G vL(v, U) or v e uL(u, U). 

(c) If v G C/ and t> G [/), then L(v, U) C L(w, [/). 

(d) If V C ?7, then L{u, U) C L(u, 1/) for all u G V. 

Elements of C/) are called multiplicative for m. If w G uL(u, U) we shall write 
and call u (L-)involutive divisor of w. The monomial w in its turn is called 
(L-) involutive multiple of u. In such an event monomial v = w/u is multiplicative for 
14 and the equality w = uv will be written as w = uxv. If w is the conventional divisor 
of w but not involutive one we shall write, as usual, w = u ■ v. Then v is said to be 
non-multiplicative for u. 

Definition 2.2 We shall say that involutive division L is globally defined if for any 
u G M its multiplicative monomials are defined irrespective of the monomial set U 3 u, 
that is, if L(u, U) = L{u). 

Definition [2.1| for every u G U C M provides the partition 

{x u ...,x n } = M L (u, U) U NM L (u, U), M L (u, U) n NMl(u, U) = (2) 

of the set of variables into two subsets: multiplicative M L (u,U) C L(u,U) and non- 
multiplicative NMl(u,U) G" L(u,U). Conversely, if for any finite set {7 G M and any 
u G Z7 the partition @ is given such that the corresponding submonoid L(u, U) of 
monomials in variables in Ml(u, U) satisfies the conditions (b)-(d), then the partition 
generates the involutive division. 

The conventional monomial division, obviously, satisfies condition (b) only in the 
univariate case. 

In what follows monomial sets are assumed to be finite, unless involutive division 
L is globally defined. In this case, since L is defined irrespective to the monomial set, 
it admits extension to infinite sets. 

Definition 2.3 A monomial set U G M is involutively autoreduced or L—autoreduced 
if the condition uL(u, U) H vL(v, U) = holds for all distinct u,v G U . 

Definition 2.4 Given an involutive division L, a monomial set U is involutive^ with 
respect to L or L— involutive if 

U ueu uM = U ueU uL(u,U). (3) 



1 Janet |l7j and Thomas [E3| call such sets complete. 
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Definition 2.5 An L— involutive monomial set U is called L— completion of a set U C 
E/ if 

U u6C/ mM = U^«L(u,C/) . 

If there exists a finite L— completion U of a finite set Z7, then the latter is finitely 
generated with respect to L. The involutive division L is noetherian if every finite set 
£7 is finitely generated. 

Proposition 2.6 ffiOj If involutive division L is noetherian, then every monomial ideal 
has a finite involutive basis U . 

Proposition 2.7 If U is a finitely generated monomial set, then so is set obtained by 
autoreduction of U in the sense of the conventional monomial division. 

Proof It follows immediately from observation that any involutive completion of U is 
also an involutive completion of its autoreduced subset. □ 

Definition 2.8 A monomial set U is called locally involutive with respect to the invo- 
lutive division L if 

(Vw G U) (Wxi G NM L (u, U)) (3v G U) [ v\ L (u ■ Xi ) } . 



Definition 2.9 A division L is called continuous if for any finite set U G M and for 
any finite sequence {iii}(i<j<fc) of elements in U such that 

(Vi < k) (3xj G NM L (u u U)) [ u i+1 \ L Ui ■ x 3 ] (4) 

the inequality ^ Uj for % ^ j holds. 

Theorem 2.10 / p?H/ If involutive division L is continuous then local involutivity of any 
monomial set U implies its involutivity. 

Definition 2.11 A continuous involutive division L is constructive if for any UcM, 
u G U, Xi G NMl(u, U) such that u ■ Xi has no involutive divisors in U and 

(Vf G U) (Wxj G NMl(v, U)) (v ■ Xj\u ■ Xi, v ■ Xj ^ u ■ Xi) [v ■ Xj G U u< zjj u L{u, U) ] 

the following condition holds: 

(Ww G U u€U uL(u, U)) [u- Xi ^wL(w,UU{w})}. (5) 
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Given a finite set of polynomials FcM and an admissible ordering multiplicative 
and non-multiplicative variables for / £ F are defined in terms of lm(f) and the leading 
monomial set lm(F). 

The concepts of involutive polynomial reduction and involutive normal form are in- 



troduced similar to their conventional analogues [[II]] with the use of involutive division 
instead of the conventional one. 

Definition 2.12 Let L be an involutive division L on M, and let F be a finite set of 
polynomials. Then we shall say: 

(i) . p is L— reducible modulo f £ F if p has a term t = au£T(a^0) such 

that u = lm(f) x v , v £ L(lm(f),lm(F)). It yields the L— reduction p — ^ g = 
p-(a/lc(f))fv. 

(ii) . p is L— reducible modulo F if there exists / £ F such that p is L— reducible 

modulo /. 

(iii) . p is in L— normal form modulo F if p is not L— reducible modulo F. 

We denote the L— normal form of p modulo F by NFi(p, F). In contrast, the conven- 
tional normal form will be denoted by NF(p,F). If monomial u is multiplicative for 
lm(f) (/ £ F) and h = fu we shall write h = f x u. 

Definition 2.13 A finite polynomial set F is L— autoreduced if the leading monomial 
set lm(F) of F is L— autoreduced and every / £ F does not contain monomials invo- 
lutively multiple of any element in lm(F). 



Theorem 2.14 j^j If set F C R is L- autoreduced, then NF L (p,F) — iff p £ 
R is presented in the form p = c ifi x u %j where fa £ F , q £ and w^- £ 
L(lm(f), lm(F)) are such that mj ^ for i ^ k. 



Corollary 2.15 j^/ If polynomial set F is L— autoreduced, then NFi(p, F) is uniquely 
defined for any p<ER, and NF L (pi +p 2l F) = NF L (p u F) + NF L (p 2 , F). 

Definition 2.16 An L— autoreduced set F is called (L-) involutive if 

(V/ £ F) (Vu £ M) [ NF L (fu, F) = ] . 

Given v £ M and an L— autoreduced set F, if there exist / £ F such that lm(f) -< i> 
and 

(V/ £ F) (Vu £ M) (Zm(/) • u -< u) [ NF L (fu, F) = ] , (6) 

then F is called partially involutive up to the monomial v with respect to the admissible 
ordering -<. F is still said to be partially involutive up to v if u -< lm(f ) for all f £ F. 
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Theorem 2.17 jffilfl An L—autoreduced set FcIm involutive with respect to a con- 
tinuous involutive division L iff the following (local) involutivity conditions hold 

(V/ G F) {Mx t G NM L (lm(f), lm(F))) { NF L (f ■ x h F) = ] . 

Correspondingly, partial involutivity (j^) holds iff 

(V/ G F) (Vx 4 G NM L (lm(f) } lm(F))) (lm(f) ■ x, t < v) { NF L (f ■x i ,F) = 0]. 



Theorem 2.18 If F CM is an L— involutive basis, then it is also a Grobner basis, 
and the equality of the conventional and L— normal forms NF(p, F) = NFi(p, F) holds 
for any polynomial If set F is partially involutive up to the monomial v, then 

the equality of the normal forms NF(p, F) = NFi(p, F) holds for any p G E such that 
lm(p) -< v. 



Theorem 2.19 f$2y Let F be a finite L—autoreduced polynomial set, and let g ■ x be 
a non-multiplicative prolongation of g G F. Then NFi(g ■ x,F) = if the following 
holds 



(yh G F) (Vu G M 
(3f,f ,g eF) 



( lm{h) ■ u -< lm(g ■ x) ) [ NF L {h ■ u, F) = ] 

lm(f )\lm(f) , lm(g Q )\lm{g) 
lm(f)\ L lm(g ■ x) , lcm(f , g ) -< lm(g ■ x) 
NF M-^ v F)=NF L (9o-§S- v F)=0 



lt(go) - 



3 Examples of Involutive Divisions 



First of all, we give three examples of involutive division used in [IB, [L7|, |23[ for analysis 
of algebraic differential equations. For the proof of validity of properties (b)-(d) in 
Definition 12.11 for these divisions we refer to IF22I . 



Example 3.1 Thomas division pcifl . Given a finite set U C I 
considered as multiplicative for u G U if degi(u) = max{degi(v) 
multiplicative, otherwise. 



the variable Xj is 
v G U}, and non- 



Example 3.2 Janet division ||17| . Let set U C M be finite. For each 1 < i < n divide 
U into groups labeled by non- negative integers d\, . . . ,df. 



[d 



1: 



cL 



{ u G U | dj = degj(u), 1 < j < i }. 



A variable Xi is multiplicative for u G U if i — 1 and deg\{u) = max{degi(v) \ v G U}, 
or if i > 1, u G [di, . . . , <ij_i] and degi(u) = max{degi(v) \ v G [di, . . . , 
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Example 3.3 Pommaret division [16]. For a monomial u = xf 1 ■ ■ • xf. k with dk > 
the variables Xj,j > are considered as multiplicative and the other variables as 
non-multiplicative. For u = 1 all the variables are multiplicative. 

Now we present two more examples of divisions which, as does Thomas division, do 
not rest on the variable ordering. 

Example 3.4 Division I. Let U be a finite monomial set. The variable x, is non- 
multiplicative for u G U if there is v G Z7 such that 

x* 1 • • • xf™w = Zcm(w, f ), 1 < m < [n/2], ci,- > (1 < j < m) , 

and x^ G {x^ , . . . , x^ m j . 

Example 3.5 Division II. For monomial u = xf 1 ■ ■ -x d k n the variable Xj is multiplica- 
tive if di = d max (u) where d max (u) = max{di, . . ., d n }. 

To distinguish the above divisions, the related subscripts T, J, P, I, II will be used. 
We note that 

• Thomas division, Divisions I and II do not depend on the ordering on the variables 
Xj. Two other divisions, as defined, are based on the ordering ([]]). 

• Pommaret division and Division I are globally defined in accordance with Defi- 



nition 2.1 , and, hence, admit extension to infinite monomial sets. 



Proposition 3.6 Divisions I and II are involutive. 



is 



Proof Division I. First of all, we prove that the condition (b) in Definition 
fulfilled. Let u ^ v be elements in U such that u\iw and v \iw for some w G M. If 
u\iv or v\ju, then we are done. Otherwise, lcm(u,v)/u or lcm(u,v)/v contains non- 
multiplicative variables for u or v, respectively. Because lcm(u,v)\w, it follows that w 
cannot be involutively multiple of both u and v. 

Consider now u G U such that u\iv for some v G U, and v u. Suppose v\iw 
for some w G M, and assume for a contradiction that w is not involutively multiple of 
u. Then there are variables Xj 1; . . . ,Xj m (1 < m < [n/2]) containing in w/v which are 
non-multiplicative for u and there is t G U such that ux^ ■ ■ ■ x\™ = lcm(u, t). Because 
v/u does not contain x ii; . . . ,Xj m it follows vx^ ■ ■ ■ x^ = lcm(v,t), that contradicts 
our assumption that w G vL(v, U) and proves the fulfillment of condition (c). 

The condition (d) holds too, since an enlargement of the set U may, obviously, only 
produces extra non-multiplicative variables for any u G U. 

Division II. Let u with d u = d max (u) be an involutive divisor of some monomial 
w G M. Then, by definition, degi(u) = min(degi(w) , d u ) (1 < i < n). Thus, given 
monomial w and number d u such that d u < d w where d w = d max (w), the corresponding 
involutive divisor u of w is uniquely defined. If there are two involutive divisors u, v of 
w with d u < d v , then it follows that 
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degi(u) = degi(v) = degi(w) if degi{w) < d u , 
d u < degiiy) = min(degi(w) , d v ) if degi(w) > d u . 

Hence, u is involutive divisor of v and the condition (b) is fulfilled. 

The condition (c) is an easy consequence of the relations deg^u) = min(degi(v), d u ) 
and degi{v) = min(degi(w) , d v ) . 

The condition (d) holds trivially, because the division as well as Pommaret one does 
not depend on monomial set U at all. □ 



Proposition 3.7 For any finite monomial set U and for any monomial u G U, the 
inclusion Mt(u,U) C Mj(u,U) and, respectively, NMj C NMt{u,U) holds. 

Proof If Xi G NMj(u,U), then, obviously, degiiu) < hi = max{degi{u) \ u G U}, 
and, hence, Xi G NMt(u,U). □ 



Example 3.8 U = {x 2 ,xy,z} (x y y y z). 



monomial 


Thomas 


Janet 


Pommaret 


Division I 


Division II 


M T 


NM T 


Mj 


NMj 


M P 


NM P 


Mi 


NMj 


M n 


NM n 


x 2 


X 


y,z 


x,y,z 




x,y,z 




X 


y,z 


X 


V,z 


xy 


y 


X, z 


y,z 


X 


V,z 


X 


y 




x,y 


z 


z 


z 


x,y 


y,z 


X 


z 


x,y 




X 


z 


x,y 



Proposition 3.9 Divisions given by Examples 3Tl-375 are continuous and construc- 



tive. All these divisions except that of Pommaret are also noetherian. 



Proof The proof for Thomas, Janet and Pommaret divisions is given in |22] . Consider 
Divisions I and II. 

Continuity. Let U be a finite set, and {ui} (i<i<M) be a sequence of elements in 
U satisfying the conditions (|j). In accordance with Definition |2.9| we shall show that 



there are no coinciding elements in the sequence for each of the two divisions. There 
are the following two alternatives: 



(0 



Ui-i ■ X i 



3 1 



III 



Ui ^ Ui-i ■ Xj 



(7) 



Extract from the sequence the subsequence = Ui k }(i<k<K<M) of those elements 
which occur in the left-hand side of relation (ii) in (ffh. 

Division I. Show that tk\ilcm(tk-i, tk) and ^ lcm(tk-i,tk). We have tk x Wk = 



u 



_i • x 



tk-i ■ Vk-i where -^Wk\vk-i- Indeed, suppose Wk\Vk-i- Apparently, this 
implies the relation tj~ = ur Z\ where i^-i < I < ik, and the variable Xj t G NMj(ui, U), 



jk 



which figures in Definition of the sequence satisfies Xj^Wk and -iXj^Zi. It 
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follows that lcm(t k , wz+i) = t k Xj t what, in accordance with definition of the division in 
Example |3.4| , contradicts multiplicativity of Xj t for t k . Therefore, we obtain the relation 

tk ■ Vk = tk+l X Wk+l , ,gN 

gcd(v k , Wk+i) = gcd(v k , w k ) = l, 

where w k+ i contains more then [n/2] variables with positive exponents, and, hence, v k 
contains only non-multiplicative variables for t k . 

We claim now that any Vj occurring in (|]) with j > k as well as v k contain only 
non- multiplicative variables for t k . For j = k + 1 we multiply t k v k by v k +i 

tkVkVk+l = (tk+l ■ V k+1 )w k+ i = (t k+2 X W k+2 )w k+1 , 

gcd(v k ,w k+1 ) = gcd(v k+1 ,w k+1 ) = gcd(v k+1 ,w k+2 ) = 1 . 
It yields 

t k V k V k+ i = (tfc+2 X W k+2 )W k+ i , , s 

gcd(v k v k+1 ,w k+2 w k+1 ) = 1. 

Because w k+ i contains more than [n/2] variables, the number of variables occurring in 
the product v k v k+ i is less or equal [n/2], and, thus, variables which are multiplicative 
for t k are not contained in v k+ \. 

If we proceed, sequentially multiplying the upper equality in (|9]) by v k+ j (j = 2, . . .), 
rewriting the right-hand side of every product in terms of t k+ j + i and cancelling the 
common factors, then we obtain the equality 

t k V k ■ ■ ■ V k+ j-\V k+ j = (tfe+j+1 X W k+ j + i)w k+ i ■ ■ ■ W k+ j-iW k+ j , 

gcd(v k ■ ■ ■ v k+j ^iv k+j , w k+ j +1 w k+1 ■ ■ ■ w k+j -iw k+j ) = 1 . 

It proves the claim and implies ti ^ tj for i ^ j. 

It remains to prove that elements of the sequence {wi}(i<i<A/) which occur in the left- 
hand side of relation (i) in (0) are also distinct. Assume for a contradiction that there 
are two elements Uj = u k with j < k. In between these elements there is, obviously, 
an element from the left-hand side of relation (ii) in (^). Let U{ m (j < i m < k) be 
the nearest such element to Uj. Considering the same non-multiplicative prolongations 
of u k as those of Uj in the initial sequence, one can construct a sequence such that 
the subsequence of the left-hand sides of relation (ii) in ([?]) has two identical elements 
u i k = u im witn i k > i m - 

Division II. The above defined elements t k which occur in the left-hand side of the 
relation (ii) in (|7|) are distinct because d ma x(t k+ i) < d ma x(t k ). The other elements 
occurring in relation (i) in (0) are also distinct since deg{ui k+ j) = deg(ui k+ j^\) + 1 
(j = 1, . . . ,i k+1 - i k - 1) and 

dmaxitk) d max (Ui k ^-i) • d max (ltj fe+1 — i) . 

Constructivity. Division I. Let u ■ x i; u e U, Xi G NMi(u,U) be a non- 
multiplicative prolongation such that 

u ■ Xi = U\v x w, U\ G U, v G I(ui, U), w G I(u\V, U U {uiv}), to / 1. 
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Show that if Xj \ w, then Xj G Mi(ui,U). Suppose a:.,- G NMi(ui,U). It means that there 
is v\ G U satisfying degj{u\) < degj{v\). Because ->Xj\v, we have degj(u%v) < degj(vi), 
and, hence, Xj G NMi(u%v, U U {u\v}). 

Division II. Since this division is globally defined, its constructivity ia an immediate 
consequence of the property (c) in Definition |2.1| . 

Noetherity. Division I. Its noetherity follows from Proposition |3.7| and noetherity 
of Thomas division, since every Thomas completion of a set U, obviously, is also its 
completion with respect to Division I. 

Division II. Given a finite set U C M and u G U with d u = d max (u), complete 
the set by the monomial xf u ■ ■ ■ xf^ and all its divisors multiple of u. If we do such a 
completion for every u G U we obtain, apparently, an involutive completion of U. □ 



4 Minimal Involutive Monomial Bases 

Let U be a finitely generated monomial set with respect to involutive division L. In this 
case a finite involutive completion U ~D U forms the involutive basis of the monomial 
ideal generated by U. A monomial ideal may not have the unique involutively autore- 
duced basis. For instance, from the definition of Janet division given in Example |3.2| it 



is easy to see that any finite monomial set is Janet autoreduced. Therefore, enlargement 
of a Janet basis by a prolongation of any its element and Janet completion of the en- 
larged set leads to another Janet basis of the same monomial ideal. Similarly, Thomas 
division and Division I do not provide uniqueness of involutively autoreduced bases 
whereas Pommaret division and Division II do, as the following proposition shows. 

Proposition 4.1 Let L be a globally defined involutive division. Then any monomial 
ideal has the unique L— involutive basis. 

Proof Assume that there are two distinct L— bases U\ and U 2 of the monomial ideal 
Id(U) where U is the finite monomial set generating the ideal and autoreduced in the 
sense of the conventional monomial division. Both U\ and U 2 are apparently involutive 
completion of U. It follows U\ \ U 2 ^ and U 2 \ U\ ^ 0. Otherwise one of sets U\, U 2 
would contain another, and, hence, could not be L— autoreduced. Indeed, let U 2 C U\. 
Then any element of u G JJ\ \ U 2 is multiple of some element in U, and, in accordance 



with Definition |27|, u is involutively multiple of some element v G U 2 . 



We obtain that for any u G U\ \ U 2 there is v G U 2 \ U\ such that v\jju and for any 



v G U 2 \Ui there is w G Ui\U 2 such that w\lv. Thus, by property (c) in Definition |27l 



given u G U\ \ U 2 there exist w G U\\U 2 such that w\lU. Since U\ is L— autoreduced, 
it is possible only if u — w. But this implies u = v. The obtained contradiction proves 
the proposition. □ 

Definition 4.2 Let L be an involutive division, and Id(U) be a monomial ideal. Then 
its L— involutive basis U will be called minimal if for any other involutive basis V of 
the same ideal the inclusion U C V holds. 
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Proposition 4.3 If U C M is a finitely generated set with respect to a constructive 
involutive division, then monomial ideal Id(U) has the minimal involutive basis. 



Proof The proof follows immediately from Proposition EO and existence of the min- 



imal involutive completion for a finitely generated set |Z| . □ 



If L is constructive, then to compute the minimal involutive basis for an ideal generated 
by a given finite monomial set one can use the following algorithm which is a slightly 



modified version of algorithm InvolutiveCompletion in paper [p2|1 . 

Algorithm MinimallnvolutiveMonomialBasis: 

Input: U, a finite monomial set 

Output: U, a minimal involutive basis of Id(U) 

begin 1 

U := Autoreduce(U) 2 

choose any admissible monomial ordering -< 3 

while exist u & U and x G NMl(u, U) s.t. 4 

u ■ x has no involutive divisors in U do 5 

choose such u, x with the lowest u ■ x w.r.t. -< 6 

U:=UU{u-x} 7 

end 8 

end 9 

The proof of correctness and termination, for a finitely generated set, of this algorithm 

is the same as that of algorithm InvolutiveCompletion p21 if Proposition 1277 



is 

taken into account. In effect, the below algorithm constructs the minimal involutive 
completion of an autoreduced, in the sense of the conventional monomial division, 
initial monomial set. This autoreduction is just done in line 2 of the algorithm. 

Example 4.4 (Continuation of Example |3l|). The minimal involutive bases of the 
ideal generated by the set U = {x 2 , xy, z} (x y y y z) are given by 

Ut = {x 2 , xy, z, xz, yz, x 2 y, xyz, x 2 z, x yz} , 

Uj = {x 2 ,xy,z,xz} , 

Up = {x 2 , xy, z, xz, yz, y 2 z, y k z, . . .} , 

Uj = {x 2 , xy, z, xz, x 2 y, xyz, x 2 z, x 2 yz} , 

U n = {x 2 ,xy,z,xz,yz,xyz} , 

where k e N (k > 2), and subscripts in the left-hand sides stand for different involutive 
divisions considered in Section 3. This example explicitly shows that Pommaret divi- 
sion is not noetherian. However, for another ordering z y x y y the set U is finitely 
generated, and then Up = U. 

One should note that selection of a L— irreducible non-multiplicative prolongation 
which is lowest with respect to an admissible monomial ordering and which we call 
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normal is of fundamental importance for the above algorithm. We demonstrate this 
fact by the following example. 

Example 4.5 Let U = {x 2 ,xz,y} and L be Pommaret division with x y y y z. By 
the normal selection strategy, the lowest irreducible non-multiplicative prolongation is 
y ■ x with respect to any admissible monomial ordering. Enlargement of U by xy gives 
the Pommaret basis U = {x 2 ,xy,xz,y} of ideal Id(U) which is obviously minimal. 
This shows that U is a finitely generated set. However, if we would take first the 
prolongation xz ■ y which is involutively irreducible modulo U, but not lowest, then we 
might obtain the infinite chain of irreducible prolongations: 

xz — > xyz — > xy 2 z ■ ■ ■ — > xy k z — > • • • 



Definition 4.6 Let L be a constructive involutive division, U be a finite monomial 
set and V = Autoreduce(U) . Then set U will be called (L-) compact if U = V or U is 
obtained from V in the course of the above algorithm. 

As an immediate consequence of this definition we have the following corollary. 

Corollary 4.7 If U C M is a finitely generated set with respect to a constructive 
involutive division L, then a compact involutive basis of ideal Id(U) is minimal. 



5 Minimal Involutive Bases of Polynomial Ideals 



In paper |22|] we proposed the next algorithm InvolutiveBasis for computation of 
involutive bases of polynomial ideals. In the algorithm the initial polynomial set F is 
subject, first of all, to the conventional autoreduction in line 2. Next are two main 
steps which are sequentially made: 

(i) . By the normal strategy, a non- multiplicative prolongation g ■ x of element g in 

the intermediate basis G with the lowest lm(g ■ x) is selected in line 5. If there 
are several different non-multiplicative prolongations with the same leading term, 
then any of them may be selected. 

(ii) . If h = NFi(g-x, G) ^ 0, then G is enlarged by h, and the involutive autoreduction 

of the enlarged set is done in line 8. 

In order to apply the criterion in line 7 for elimination of superfluous involutive reduc- 
tions and also to avoid repeated prolongations, the auxiliary set T of triples (g, u, P) 
is used. Here g G G, and u is either the lowest, with respect to the ordering -<, leading 
monomial in lm{G) such that g was produced by non-multiplicative prolongations of 
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/ G G with u = lm(f), or u = lm(g) if there is no such / in G. Those variables in 
NMi(g, G) have been chosen in line 5 are collected in set P. 

Algorithm InvolutiveBasis: 

Input: F, a finite polynomial set 

Output: G, an involutive basis of the ideal Id(F) 

begin 1 

G := Autor educe(F); T : = 2 

for each g G G do T := T U {(g, lm(g), 0)} 3 

while exist (g, u, P) G T and x G NM L (lm(g), lm{G)) \ P do 4 

choose such (g,u,P),x with the lowest lm(g) ■ x w.r.t. -< 5 

T:=T\{( fl)1 i,P)}U{( ft! i,PUW)} 6 

if Criterion(g ■ x, u, T) is false then h := NF^(g ■ x, G) 7 

if h ^ then G := Autor educe L (G U {/i}) 8 

if lm(h) = lm(g ■ x) then T:=TU {(/i, w, 0)} 9 

else T := T U {(/i, lm(h), 0)} 10 

Q:=T; T := 11 

for each g G G do 12 

if exist (f,u,P) G Q s.t. lm(f) = lm(g) then 13 

choose gi E G s.t. Zm(gi)|iW 14 

T:=TU{(g,lm( gi ),P)} 15 

else T:=TU{(j, lm(g), 0)} 16 

end 17 

end 18 

Criterion(g, u, T) is true provided that if there is (/, v,D) G T such that lm(f)\Llm(g) 
and lcm(u,v) -< lm(g). Correctness of this criterion, which is just the involutive 



form p2| of the Buchberger's chain criterion, is provided by Theorem 2.19 



Definition 5.1 Given a constructive division L, a finite involutive basis G of ideal 
Id(G) is called minimal if lt(G) is the minimal involutive basis of the monomial ideal 
generated by {/*(/) | / G Id(G)}. 



Theorem 5.2 A monic minimal involutive basis is unique. 

Proof Assume for a contradiction that a polynomial ideal Id(F) has two distinct monic 
minimal involutive bases G\ and G2- Their minimality means that lm{Gi) = lm{G2). 
Since Gi and G 2 are distinct there are g\ G G\ and g 2 G G 2 such that lt(gi) = U{g2) 
but gi 7^ g 2 - Since g\ — g 2 G Id(F), by Theorem |2.18| , we have NF L (gi — g 2) G\) = 
NF L (gi — g 2 ,G 2 ) = 0. Therefore, at least one of the sets Gi,G 2 is not involutively 
autoreduced, and, hence, in accordance with Definition |2.16| , it cannot be involutive 
basis. □ 
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For a globally defined involutive division, by Proposition |4.1| , this proof, obviously, 
is also valid for polynomial ideals with infinite involutive bases. Therefore, we have the 
following corollary. 



Corollary 5.3 Given a globally defined involutive division, every polynomial ideal has 
the unique involutive basis. 

Thus, given a globally defined involutive division L, the output of algorithm Invo- 
lutiveBasis, in the case of its termination, is unique for a given polynomial ideal 
irrespective of an ideal generating set F in the input. 

However, even though the algorithm may not terminate it is still able to compute 
a Grobner basis as the following proposition shows. 

Proposition 5.4 Let L be a continuous involutive division and G be an intermediate 
polynomial basis generated by algorithm InvolutiveBasis. If the ordering -< is degree 
compatible, then in a finite number of steps G becomes a Grobner basis. 

Proof Let the current prolongation g ■ x is such that h = NFi(g ■ x,G) ^ 0. Then at 
the second main step of the algorithm (step (ii) as described above), the intermediate 
polynomial set is enlarged by h. In so doing there are two alternatives: 

(a) lm(h) = lm(g ■ x); (b) lm(h) -< lm(g ■ x) . 

In the latter case lm(g ■ x) is involutively reducible by some lt(f) G lt(G), that is, 
lm(g) ■ x = lm(f) x w. Then, by Theorem |2.14| and Corollary |2.15| we have the equality 
NF L (g ■ x, G) — NF L (S(f, g), G) where S(f, g) = g ■ x — f x w is an S— polynomial. 

In this case, unlike the case (a), the monomial ideal Id{lm(G)) is changed. Indeed, 
let there is a polynomial h\ G G such that lm(h) is multiple of lm(h\) but not invo- 
lutively multiple, that is, lm(h) = lm{h\) ■ (lm(h) / 'lm(hi)) . By the normal selection 
strategy, set G satisfies the condition (|6]) of partial involutivity up to the monomial 
lm(h) with respect to the ordering -< what implies NF^h, F) = 0. 

Furthermore, by Theorem JFRj NF L (S(g u g 2 ),G) = NF(S(g u g 2 ),G) = for any 
S— polynomial S(gi,g2), {g\-,g% G G) with lcm(lm(gi) , lm(g2}) -< lm(g ■ x). 

It remains to prove that every S(gi,g2) such that NF(S(gi, g 2 ),G) ^ is com- 
puted at some step of the algorithm. Since set G is L— autoreduced, monomial u = 
lcm(lm(gi) , lm(g 2 )) cannot be involutively multiple of both lm(gi),lm(g 2 ). Hence, by 
degree compatibility of the ordering -<, in a finite number of steps at least one of 
g±, g 2 will be non-multiplicatively prolonged to a polynomial g with lm(g) = u. Let 
g be obtained by non-multiplicative prolongations of g±, and the current prolongation 
be g with u = lm(gi) ■ (u/lm(gi)). If u is involutively multiple of lm(g 2 ) or Im(g^) 
where g^ is a polynomial obtained in the course of the algorithm by non-multiplicative 
prolongations of g 2 , then we are done. 

Otherwise, there is to be g G G such that u = lm(g) = lm(g 2 ) ■ (u/lm(g 2 )), and one 
of the two polynomials g,g will be constructed before the another. Since their leading 
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monomials coincide, the leading monomial of the latter will be involutively reducible 
by the leading monomial of the former. □ 

Though, by Corollary |5.3| , algorithm InvolutiveBasis, if it terminates, computes 
the minimal involutive basis for a globally defined involutive division it may not be 
the case for arbitrary involutive division. If we use, for instance, any of divisions in 
Examples |3.1[ - |3.2| and |3.4| , then, given a polynomial ideal Id(F), the algorithm output 
depends on the structure of input generating set F. 

Example 5.5 Let F = {x 2 y — 1, xy 2 — 1, y 4 — 1}. The lexicographical Janet basis for 
x y y y z computed by algorithm InvolutiveBasis is 

{ x 2 y - 1, x 2 - 1, xy 2 - 1, xy - 1, x - 1, y* - 1, y 3 - 1, y 2 - 1, y - 1 } . 

The reduced Grobner basis {x — l,y — 1} of Id(F) is also the minimal Janet basis. 



Proposition 5.6 // algorithm InvolutiveBasis takes a reduced Grobner basis as in- 
put it produces a minimal involutive basis for a constructive involutive division. 

Proof Let g-xbe & non-multiplicative prolongation of element g in intermediate poly- 
nomial set G, and h = NFi(g ■ x,G). We note that either h = or lm{h) = lm(g ■ x). 
Otherwise, as shown in the proof of Proposition |5.4| , lm(h) would not belong to mono- 
mial ideal Id{lm(G)) = Id(lm(F)). Thus, the output monomial set lm{G) is con- 
structed just as it would be done by applying algorithm MinimallnvolutiveMono- 
mialBasis to lm(F). It follows that lm(G) is the minimal basis of Id(lm(F)). □ 

The next algorithm constructs a minimal involutive basis, and generally deals with 
less number of intermediate polynomials than algorithm InvolutiveBasis causing the 
computational efficiency to increase. 

Theorem 5.7 Let F be a finite subset o/R and L be a constructive involutive division. 
Suppose ordering y is degree compatible. Then algorithm MinimallnvolutiveBasis 

computes a minimal involutive basis of Id(F) if this basis is finite. If L is noetherian, 
then the basis is computed for any ordering. 

Proof Correctness. First of all, we recall that correctness of the involutive criterion 



which is verified in lines 14, 23 follows from Theorem 2.19 . As distinct from the 
algorithm InvolutiveBasis here are two disjoint subsets T and Q of the triples. They 
are built in such a way that lm(g) -< lm(f) for any g in (g,u,P) G T and / in 
(f,v,D) G Q. Let G be a polynomial set {g \ (g,u,P) G Q}. First of all, we claim 
that ideal Id(G U G) is an invariant of the repeat-loop. Indeed, it is trivially true 
upon initialization. Inside the loop, if a polynomial is removed from G in lines 18 and 
28, then it is added to G. On the other hand, removal of a triple from Q, that is, the 
corresponding polynomial from G in line 11, does not change G iff NFi(g, G) = 0. 
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Furthermore, set T is handled by the lower while-loop in lines 19-29 just as it done 
in algorithm InvolutiveBasis except for restriction in line 20 and the set contraction 
in lines 27-28. In the latter case all the elements in G with lm(g) y lm(h), where h 
is the normal form of the current prolongation, are moved to G. Thus, this while- 
loop preserves the property of partial involutivity up to monomial v -< lm(h) for the 
intermediate set G, in accordance with Theorems [2.17| and [2.18] , if there is a partially 



involutive set in the input of the loop. Besides, two elements with coinciding leading 
terms obviously never occur in set G. 

Algorithm MinimallnvolutiveBasis: 

Input: F, a finite polynomial set 

Output: G, the minimal involutive basis of the ideal Id(F) 

begin 1 

F := Autoreduce(F) 2 

choose g G F with the lowest lm(g) w.r.t. -< 3 

T:={(g,lm(g)M; Q 0; G := {g} 4 

for each / G F\ {g} do 5 

Q:=QU{(f,lm(f)M 6 

repeat 7 

h:=0 8 

while Q ^ and h = do 9 

choose g in (g,u,P) G Q with the lowest lm(g) w.r.t. -< 10 

Q:=Q\{(g,u,P)} * ' 11 

if Criterion(g,u,T) is false then h:=NF L (g,G) 12 

end 13 

if h ^ then G := G U {h} 14 

if lm{h) = lm{g) then T := T U {{h, u, P)} 15 

else T:=TU {(h, lm(h), 0)} 16 

for each / in (f,v,D) G T s.t. lm(f) y lm(h) do 17 

T:=T\ {(/,«,£>)}; Q := Q U {(/, «, D)}; G :=G \ {/} 18 

while exist (c/, u, P) G T and ac G NM L (g, G)\P and, if Q ^ 0, 19 

s.t. Zm(p ■ x) -<; /m(/) for all / in (/, v,D) G Q do 20 

choose such (g,u,P),x with the lowest lm(g) ■ x w.r.t. -< 21 

T:=T\{(g,u,P)}U{(g,u,PU{x})} 22 

if Criterion(g ■ x, u, T) is false then h := NFi(g ■ x, G) 23 

if /i^O then G:=GU{/i} 24 

if lm{h) = lm{g ■ x) then T := T U {(/i, u, 0)} 25 

else T := T U {(/i, lm(h), 0)} 26 

for each / in (/, v,D)eT with lm(f) y lm(h) do 27 

!T:=T\ {(/,«,£>)}; Q := Q U {(/, «, D}); G := G \ {/} 28 

end 29 

until Q^0 30 

end 31 
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In what follows polynomials in G, if G ^ 0, are successively selected in accordance with 
the normal strategy; taken out of the set and L— reduced modulo G. The upper while- 
loop in lines 9-13 proceeds until the normal form h of the selected polynomial does not 
vanish. Then set G is enlarged by h in line 14. The repeat-loop terminates when set 
G becomes empty in line 11 and the lower while-loop does not lead to appearance of 
new elements in this set. It means that the output set G is an involutive basis of ideal 
Id(G) = Id(F). 

Now, by Corollary 3~7 , to prove minimality of the output basis it is sufficient to show 
that the lower while-loop always ends up with L— autoreduced polynomial set G such 
that lt(G) is compact. As we have already seen, this loop preserves partial involutivity. 
Initially there is a single polynomial which has the minimal leading monomial, and, 
therefore, its handling in the loop produces a compact leading monomial set. 

Suppose a partially involutive polynomial set G with compact lm(G) was produced 
by the lower while-loop, and then it is enlarged by h = NFi(g, G) in line 14 when G 
is partially involutive up to some monomial v -< lm(g). 

If lm(h) = lm(g), then, by restriction in line 20, lm(h) >- lm(f) for all / G G. 
By property (d) in Definition [2.1| , we obtain that NMi(lm(f),lm(Gi = G U {h})) C 
NM L (lm(f), lm(G)) for any f £ G. 

Let lm(h) has no conventional divisors in lm(G). Then, starting with the set Go — 
Autoreduce(Gi), and completing Go with irreducible non-multiplicative prolongations 
of its elements by the normal strategy, we construct set G2 ^ G± partially involutive 
up to the monomial v and with compact ImiGi). If we start now with set lm(Gi) 
and complete it, if necessary, with irreducible non- multiplicative prolongations of its 
elements in order to obtain a partially involutive set up to v, then we arrive at the same 
set G2. Indeed, even in the presence of extra intermediate elements, if G%\Gi 7^ 0, there 
cannot occur reduction of an element p G G either by an element in G or by an extra 
element. The former reduction is impossible by property (d) of involutive division. The 
latter reduction, if it would hold, by properties (c)-(d) and by Theorem |2.18| , would 
lead to reducibility of p in the earlier set G when h has not been added yet. 

If lm(h) is multiple of some element in lm(G), then continuation of processing with 
G\ in the lower while-loop yields a partially involutive polynomial set up to lm(h). 
In doing so, h is involutively reduced either to zero, or to a polynomial which changes 



the monomial ideal Id(lm(G)), as we have shown in the proof of Proposition [574 



Correspondingly, G, after contraction in lines 27-28, is reset to the partially involutive 
form with the compact leading monomial set. 

In the case when lm(h) -< lm(g), the elimination which is done in line 18 converts, 
apparently, the situation into one of two alternatives we have just considered. 

Thus, the repeat-loop, if it terminates, ends up with an involutive set G with 
compact lm(G), that is, with the minimal involutive basis. 

Termination. As it shown in the proof of Proposition 5A, there may be a finite 
number of cases when polynomial g chosen in lines 10 or prolongation g-x chosen in line 
21 have reducible leading monomials. It implies finitely many redistributions of triples 
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between T and Q done in lines 18 and 28. If Id(F) has the finite minimal involutive 
basis, and ordering -< is degree compatible, then the lower while-loop terminates 
irrespective of Q is empty set or not. This follows immediately from Propositions |5.4|, 



|5.6| and compactness of lm{G). Since the upper while-loop is obviously terminates, 
and set Q is refreshed finitely many times, in a finite number of steps the algorithm 
arrives at Q = in line 30. 

If involutive division L is noetherian then the algorithm terminates for any ordering 
-< because the lower while-loop terminates for the same reason as the while-loop does 
in algorithm InvolutiveBasis EI| . □ 



6 Conclusion 

As we noted above, algorithm MinimallnvolutiveBasis deals, generally, with less 
number of intermediate polynomials then algorithm InvolutiveBasis. Besides, if invo- 
lutive division L is not globally defined, then we may not obtain the minimal involutive 
basis in the output of the latter algorithm. But even for globally defined divisions the 
former algorithm avoids the involutive autoreduction done in the latter algorithm at ev- 
ery step of the intermediate set enlargement. That is why we expect higher efficiency 
of algorithm MinimallnvolutiveBasis with respect to algorithm InvolutiveBasis 
for arbitrary involutive division. 

One could also construct the minimal involutive basis by computing the reduced 
Grobner basis and then enlarging it by non-multiplicative prolongations of its ele- 
ments until the leading monomial set becomes involutive. To construct the reduced 
Grobner basis one can use the Buchberger algorithm or perform the conventional au- 
toreduction of an involutive basis computed by algorithm InvolutiveBasis. However, 
unlike Buchberger algorithm, algorithm MinimallnvolutiveBasis benefits from the 
involutive technique, and as we have argued is favored over the use of algorithm Invo- 
lutiveBasis for intermediate computation. 

In paper EJI] for constructing Janet bases for linear partial differential equations one 



more algorithm is described which goes back to the original computational scheme |T7 . 
Its analog in commutative algebra contains two basic subalgorithms which are succes- 
sively performed: completion of a polynomial set by non-multiplicative prolongations 
of its elements until the set of leading monomials becomes involutive or complete (see 
footnote at page 4); the conventional autoreduction of the obtained set. In this case 
due to the second subalgorithm the output Janet bases are minimal. However, such an 
algorithmic procedure is far short of optimum from the computational point of view. In 
so doing one has to perform the repeated prolongations and deal with all the possible 
S— polynomials. In our algorithm MinimallnvolutiveBasis the repeated prolonga- 
tions are eliminated by storing in the triple sets T and Q those non-multiplicative 
variables which have been used for a given polynomial. Furthermore, the use of the 
involutive analogue of the Buchberger's chain criterion allows one to cut considerably 
the number of computed S— polynomials. 
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The algorithms described in this paper just as Zharkov and Blinkov algorithm can 

and also 



be extended to systems of linear systems of partial differential equations [25 



to some classes of nonlinear systems. Being uniquely defined, minimal involutive bases 
much like reduced Grobner bases can be considered as canonical ones for polynomial 
and differential ideals. The corresponding form of partial differential equation systems 
is just the standard one. By transforming a given system into this form one can 



determine the dimension of the solution space and a set of initial conditions providing 
the existence of a uniquely defined and locally holomorphic solution [[H], |26], ^7], |28fl . 
Involutive algorithmic ideas may be also rather fruitful in constructing the canonical 



bases for finitely generated ideals in free Lie algebras and superalgebras [2S 
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